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Reprocessed emission from warped accretion discs induced 
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ABSTRACT 

The broad Balmer emission-line profiles resulting from the reprocessing of UV/X-ray 
radiation from a warped accretion disc induced by the Bardeen-Petterson effect arc 
studied. We adopt a thin warped disc geometry and a central ring-like illuminating 
source in our model. We compute the steady-state shape of the warped disc numeri- 
cally, and then use it in the calculation of the line profile. We find that, from the outer 
radius to the inner radius of the disc, the warp is twisted by an angle ~ 7r before being 
flattened efficiently into the equatorial plane. The profiles obtained depend weakly on 
the illuminating source radius in the range from 3r g to 10r g , but depend strongly on 
this radius when it approaches the marginally stable orbit of an extreme Kerr black 
hole. Double- or triplet-peaked line profiles are present in most cases when the illumi- 
nating source radius is low. The triplet-peaked line profiles observed from the Sloan 
Digital Sky Survey may be a "signature" of a warped disc. 
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1 INTRODUCTION 

A small fraction of active galactic nuclei (AGNs) with broad 
lines show double-peaked emiss ion- line profiles, which sug- 
gests an accretion di s c orig in dEracleous fc Halpern 19941. 
2003 : IStrateva etafl 12003 : iGezari. Halpern fc Eracleousl 
20071 ). In a stationary circular relativistic disc model, 



Doppler boosting will make the blue peak of the profile 
higher than the red one. However, observations of the 
variability of line emission in some objects show that at 
least at some epochs the profile asymmetry is reversed, 
with the red peak higher than the blue one, contrary 
to the predictio ns of homogeneous, cir c ular relativis- 
tic d is c models dMiller fc Peterson! Il990l : IStrateva et al.l 
120031 : IGezari et alj 120071 '). To adjust the theoretical 
to the observed profiles, non-axisymmetric or warped 
discs are usually required. Various physically plausible 
processes in the accretion flow have been considered, 
such as spiral shocks, eccentric di scs, bipolar outflows 



or a binary black h ole syst e m (IChakrabarti fc Wiita 



1994 : Eracleous et al. 19951: Storchi-Bergmann et al 



19971 : IZheng. Sulentic fc Binettel Il990l : I Veilleux fc Zheng 



199ll: iBegelman. Blandford fc Reesl Il980l ; iGaskelll Il996l : 
Zhang. Dultzin-Hacvan fc Wand l2007blh as well as a 
warped accretion disc ( Wu. Wang fc Donel l2008h . This 
warped disc makes it possible for the radiation from 
the inner disc to reach the outer parts, enhancing the 
reprocessing emission lines. Thus, warping can both provide 
the asymmetry required for the variations of the emission 
lines and naturally solve the long-standing energy-budget 
problem, because the subtending angle of the outer disc 
portion to the inner one is increased by the warp. 

There are strong observational and theoretical grounds 
for believing that accretion discs around black holes 
may be warped and twisted. Evidence for the ex- 
istence of warped discs in astr ophysical systems has 
been found from observations ( Hcrrnstc in et al. I 120051: 
Caproni. Abraham fc Mosquera Cu csta 2006; Caproni et al.l 
20071 : iMartin. Tout fc Pringld l200ct ). From a theoreti- 
cal point of view, four main mechanisms for excit- 
ing/maintaining warping in accretion discs have been pro- 
posed, namely warping th at is tidally induced by a com- 
panion in a binary s ystem (|Terquem fc Bertoutlll993l . 1 19961 : 
lLarwood etallll996h. r adiation-driven or self-induced warp 
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ing l|Pringlelll996l.ll997l:lMalonev. Begelman fc Pringlelll996l; 



Malonev fc Begelman! [l997l : iMalonev. Begelman fc Nowakl 
199ct) . magnetically driven warping 



LaH 1 19991 , 120031 : 



Pfeiffer fc Laj 120041 1. or warping driven by frame drag- 
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ging dBardeen fc Pettersonl Il975l ; iKumar fc Pringld 1 19851 ; 



lArmitage fc Nataraian I Il999f ) Note 



also that warps 
generated by gravitational interactions have been in- 
vestigated in th e liter ature. In the galactic context, 
iHunter fc Toomrel (|l969h studied the linear bending waves 
of a s elf-gravitating, isolated, thi n disc . On nuclear disc 
scales, IPapaloizou Terquem fc Linl (1 19981 ) studied the evo- 
lution of a thin self-gravitating viscous disc interacting 
with a massive object orbiting the centr al black hole. 
lUlubav-Siddiki. Gerhard fc Arnaboldil (|2009h showed that 
highly warped discs near black holes can persist for a long 
time without any persistent forcing other than by their self- 
gravity. As a result, the accretion disc in some AGNs may 
be non-planar. 

To solve the energy budget for the emission-line region 
for the disc model, reprocessing of the ultraviolet(UV)/X- 
ray continuum from the inner accretion disc is required. 
However, only a very small fraction of radiation from the 
inner disc is expected to intercept the outer part of the 
disc in the case of a fiat geometrically thin disc. Two 
types of processes have been proposed to increase the 
fraction of light incident on the d isc emission-line region. 
IChen. Halpern fc Filippenkd (|l98Sh proposed that the in- 
ner part of the accretion disc is hot and becomes geo- 
metrically thick as a result of insufficient radiative cool- 
ing, and the X-ray emission from such an ion-supported 
torus is responsible for such energy input. However, double- 
peak emitters do not always have a l ow accretion rat e 
jZhang. Dultzin-Hacvan fc Wand [2007al ; iBian et al.ll2007h . 
In contrast. ICao fc Wand [ 20Qt ) proposed that a slow- 
moving jet can scatter a substantial fraction of UV/X-ray 
light from the inner accretion disc back to the outer part of 
the disc. Although this may be likely for radio-loud double- 
peaked emitters, the majority of double-peaked emitters are 
radio-quiet. In this work, we consider a system with an ac- 
cretion disc around a Kerr black hole. The black hole spin 
is misaligned with the outer parts of the disc, as there is 
no physical reason to suppose that the angular momentum 
of the accreting mass and the angular momentum of the 
spinning black hole are always aligned. Frame dragging pro- 
duced by a Kerr black hole causes the precession of the orbit 
of a particle if its orbital plane is inclined with respect to 
the equatorial plane of the black hole. This is known as the 
Lense-Thirring effect. The combined action of the Lense- 
Thirring effect, which tends to twist up the disc, and the in- 
ternal viscosity of the accretion disc, which tends to smooth 
it out, forces the alignment between the angular momenta 
of the Kerr black hole and the accretion disc. This is known 
as the Bardeen-Petterson effect. This effect tends to affect 
the inner regions of the disc owing to the short range of 
the Lense-Thirring effect(cx r -3 ), whereas the outer part of 
the disc tends to keep its original inclination. The transition 
radius between these two regimes is known as the Bardeen- 
Petterson radius Rbp (or warp radius R w ), and its exact 
location depen ds mainly on the physical properties of the 
accretion disc 



i ds mainly on the physical properties ol the 
dBardeen fc Pettersonl 19751; IScheuer fc Feiierl 



1 19961 ; iNelson fc Papaloizoul |2000| ; IFragile fc Anninosil2005h 

Basically, the Bardeen-Petterson radius is determined by 
comparing the time-scale related to the Lense-Thirring pre- 
cession with that of warp propagation through the disc. 
We consider here the propagation of warps in thin Kep- 
lerian discs in the case for which the disc is sufficiently 



viscous that the warp propagates in a diffusive manner. 
The twisted disc itself is treated using a non-relativistic ap- 
proach with an additional term describing the gravitomag- 
netic precession of the disc rings. Both analytic and numer- 
ical steady-state solutions have shown that the evolution 
of a misaligned disc arising from the Bardeen-Petterson ef- 
fect usually produces an inner flat disc and a warped tran- 
sition r egion with a smooth gradient in the tilt and twist 
angles (|Scheuer fc Fc ilcr 1996;; lNataraian fc Armitagelll999l ; 
iMartin Pringle fc Tout! 120071 ) . If the disc is thick and/or 
its viscosity is low, the warp propagates in a wave-like 
rather than a diffusive manner (|Nelson fc Papaloizoul I2OO0I ; 
iLubow. Qgilvie fc Pringlel I2002T I. In this case, the twisted 
disc is not necessarily aligned with the equatorial plane 
at small scales. The inclination angle of a low-viscosity 
disc may even oscillate or counteralign close to the black 
hole (iLubqw et al.ll20o3 : iKing et alj|2005l ; lLodato fc Pringld 
120061 l2007l ~ 

The effect of a Bardeen-Petterson disc on iron line 
profiles treated relativistica ll y ha s been examined by 
IFragile, Miller fc Vandernootl (120051 ). and iBachevi (|l999T I 
used a non-relativistic treatment to study the broad-line 
H/3 profiles from a warped disc. A relativistic treatment of 
the effect of a warped disc with a parametrized geometry 
on the Balmer lines that can be comp ared with observa- 
tions was presented bv lWu et al] |2008). Following on from 
these calculations, the disc warping induced by the Bardeen- 
Petterson effect and its influence on the broad Balmer emis- 
sion owing to the reprocessing of the central high-energy 
radiation is investigated. In Section[3]we summarize the as- 
sumptions behind our model and present the basic equations 
relevant to our problem. We present our results in Section |3j 
and the conclusions and a discussion in Section U 



2 ASSUMPTIONS AND METHOD OF 
CALCULATION 

In this paper, we focus on how a warped disc resulting 
from the Bardeen-Petterson effect affects optical emission- 
line profiles, such as Ha and H/3. A geometrically thin disc 
and a central ring-like illuminating source (radius denoted 
by R r ) around a black hole are assumed. A short review of 
the basic equations releva nt to our proble m is given below, 
and readers are referred to lWu et alj (|2008l ) for a description 
of our numerical scheme in more detail. 



2.1 Geometrical considerations 

The disc can be treated as being composed of a series of 
concentric rings of width dR and mass 2nYiRAR at radius R 
from the central point mass M with surface density T.(R,t) 
at time t and with angular momentum L = (GMi?) 1 ^ 2 EZ = 
LI, and lying in different planes. The rings interact with 
each other through viscous stresses. Each ring is defined by 
two Eulerian angles, /3{R,t) and ^{R,i), at radius R. At 
each radius R from the center, the disc has a unit tilt vector 
l(R,t) th at varies conti nuously with radius R and time t. 
Following IPring lc (1996), the vector l(R,t) is given by 

I = (cos 7 sin /3, sin 7 sin /3, cos 0). (1) 

We define the normalized vector towards the observer as 
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iobs = (sini, 0, cos i), 



(2) 



where i is the angle between the line of sight, and the normal 
to the equatorial plane lies in the XZ plane. 

We define the coordinates on the surface of the disc as 
(R, </>) with respect to a fixed Cartesian coordinate system 
(x,y,z), where (f> is the azimuthal angle measured on the 
disc surface in the direction of flow, with (f> — n/2 at the 
ascending nodeQ The element of surface area is 

dS = [l + {R0 cos <f> + Rj' sin j3 sin <j>) e R ] R dR d<j>, (3) 

where the primes indicate differentiation with respect to R, 
and e_R is the radial unit vector. The element of radiation 
flux dF received from the ring source by unit area is 

dS- 



dF cx 



IdS 



P Pr/ r (V r )d£3r 



(4) 



where Ir(vr) is the specific intensity measured by an observer 
corotating with the ring source, p is the 3-momentum of an 
incident photon from the source, g r is a factor to describe 
the shift of photon frequency along its path, which is equal 
to the ratio of the observed frequency from the illuminated 
disc to the emitted frequency from the ring source(see be- 
low) . dSl r is the element of the solid angle subtended by the 
image of the ring source observed from the illuminated disc. 
The Balmer lines in the illuminated surface are formed as a 
result of photoionization by UV/X-ray radiation from a ring- 
like source. For simplicity, we assume that the Balmer line 
intensity is proportional to t he radiation being interc epted 
by the disc (see e.g. Fig. 6a in lCollin fc Dumondlll989l ), and 
thus the line emissivity e on the disc surface can take the 
form eoc J dF. This is an approximate assumption. A de- 
tailed treatment of line emission requires solving the vertical 
structure as well as the radiative transfer in the disc explic- 
itly, which is beyond the scope of this paper. We also assume 
that the line emission is isotropic in the comoving frame. 

2.2 The basic equations for a warped disc 

The dynamics of warping accretion discs have been dis- 
cusse d by a number of au t hors jP apaloizou "fc Pringle| 
19831; iPapaloizou" fc Lin! Il995l ; iDemianski fc Ivanovl 1 19971 ; 
Navakshinll2005l ), and the static low-viscosi t y con figurations 
were first calculated bv llvanov fc Illarionovl (|l997T ). This sec- 
tion contains a short summary of the equations describing 
the structure of the accretion disc models consid ered in this 
work. We follow the formalism of iPringlel (|l992T ) but add a 
term to describe the Lense-Thirring precession, to give 

E = -%(REVk)' 

L = -^(RVnLY+^T^' + flpxL (5) 
T vis = R'^^fl'l + q-m') 

(|Chen. Wu fc Yuan! [20091) . Here we use a dot to denote 
d/dt, and the prime symbol (') to denote d/dR. There are 
two viscosities, v\ and U2, where v\ corresponds to the az- 
imuthal shear, the standard shear viscosity in a flat disc, 
and U2 is the viscosity associated with vertical shear mo- 
tions, which smoothes out the war ping. The Lense-Thirrin, 
precession fi p is given by (see e.g. iKumar fc Pringldll98^ 
i~ip = cjp/R 3 , with: 

1 This definition of <f> differs by tt/2 from that of lPrinelj il996l) . 



2GJ , „ (GM\ . 

— and J = acM^—jj, 



(6) 



where J = Jj is the angular momentum of the black hole. 

We take both viscosities to have the same power law 
form, so that 



v\ = f 10 



(-)' 

\RoJ 



and 



( R \" 



(7) 



where i^io, ^20 an d rj axe all constants and Ro is some fixed 
radius. Although the selection of _Ro is in principle arbitrary, 
it was chosen to be the Bardeen-Petterson radius Rbp(R v ), 
which is defined as u; p /i/2 in this work. For a steady-state 
warped disc, the disc shape depends only on the radius 7? w . 
The exact value of J ? w is model-dependent and i s subject to 
some uncertainties (|Nataraian fc Prini del 1 19981 ; iKing et al.l 



some uncertainties (INataraian &£ rrmg 
120051 ; IVolonteri. Sikora fc Lasotal 120071) . 



_ 200 7J). To determ ine the 
value of -R w , either the IShakura fc Sunvaevl (1 19731 ) or the 
ICollin fc Dumondl (|l990|^ model is used. Adopt ing the disc 
properties derived by Collin fc Dumondl (|l990l ). the warp 
radius in te rms of the Schwa rzschild radius of the hole, R a , 
is given by (|King et al.ll2005l ) 



~R~ S 



6 N 1 / 4 

01 or 



O.ILe 



-1/4 



M, 



1/8 



VO.03/ 



1/8 



(—) 

V0.03/ 



-5/8 



5/8 



(8) 



Here, e is the efficiency of the accretion process (i.e. L = 
eMc 2 ), L is the accretion luminosity, Lb is the Eddington 
limit, Ms is the mass of the black hole in units of 10 s Mq, 
and a is the (dimensionless) spin of the black hole. 

Assuming a Shakura-Sunyaev d isc ("middle region" ), 
the warp radius can be expressed as (|Volonteri et al.ll2007l ): 



Rw 
~Rs 



= 3.6 x 10V /8 M„ 1/S 



M<? 



-1/4 



-5/8 



-1/2 



(9) 



or, adopting the vis cosity coefficient of a Shaku ra-Sunyaev 
disc as presented by iFrank. King fc Raind (|2002l ) , we have 



-/?„ 

IT, 



1.41 x 10 



,ao3 



-16/35 / Mc - 



-6/35 



xM 



4/35 



(10) 



For the typically expected standard values of AGNs, 7? w is 
hundreds to thousands of R s . 



2.3 Photon motions in the background metric 

We review properties of the Kerr metric and formulae for 
its particle orbits, and summarize here the basic equations 
relevant to this paper. I n Boyer-Lindquist co ordinates, the 
Kerr metric is given by (|Chandrasekharj|l983l ): 



_e^ dt 2 + e^(d0-^dt) 2 + |dr 2 



ds 2 = 
where 

e 2 " = EA/A, e 2 * = sin 2 0A/E, u = 2Mar/A, 
E = r 2 + a 2 cos 2 6, A = r 2 + a 2 - 2Mr, 
A= {r 2 + a 2 ) 2 



me 2 , (11) 
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Here M, a are the black hole mass and specific angular mo- 
mentum, respectively. The equation of motion governing t he 
orbital trajectory in the r#-plane is (|Bardeen et al.lll972l ) 



dr 



AO 



(12) 



where r e and 8 e are the starting values of r and 6. The 
^-coordinate along t he trajectory is calculated by ()Wilkinsl 
ll972l ; IViergut3ll993f ) 



d(j) 



0a S111 



+ 



r (2ar - \a 2 )dr 



(13) 



For photons emitted from the broad line region with 
radii >100r s propagating to infinity, we can neglect terms 
of order 1/r 2 and higher in equations (|12|l and ()13|1 . the 
integral over 9 can be worked out in terms of a trigonometric 
integral 



AO 



/2 



Adtf 



tv/2 sin 



±sin 



_1 (wW), 



^(1-M 2 ) 



(14) 



(15) 



Where fj, = cos 8 and ^ fi < fi + , p+ = \/q 2 /{X 2 + g 2 ). 
The integral over r can be wor ked out using inverse Ja- 
cobian elliptic integrals (see e.g. ICadez. Fanton fc Calvarul 
1 19981 ; IWu fe Wan j|2007l ) . 

The specific flux d ensity F {v ) at frequency u a is given 
by (|Cunninghamlll975l ) 



F {v ) 



g I e (v e )AQ obs 



. eg 4 5(v - ff^ c ) ^ A ' g | dr Ag. (16) 

where df2 b s is the element of the solid angle subtended by 
the image of the disc on the observer's sky and g is a factor 
to describe the shift of photon frequency along its path and 
is equal to the ratio of the observed to the emitted frequency, 
the emissivity e in the integrand is calculated by 



— |-^j-s r / r (i/ r )dS2 r 



(17) 



where g T is the redshift for radiation between the emitter's 
and the observer's frame in the special c ase of two orbiting 
systems, and is given by (IViergutd ll99ah 



9t 



7ded"(l - £2 d A) 



(18) 



7 r e r "(l - fl r X) ' 

where "f a \and'y T are the Lorentz factor measured in the locally 
non-ratating frame(LNRF). The element of solid angle A£l T 



(19) 



r d r l d(X,q) d(r,g r ) 

where is the distance from the illuminated point on the 
disc to the black hole . The two impact param e ters a and 
/3, first introduced by ICunningham fc Bardeenl (|l973l ). are 
defined as 



rp 



if) 



and (3 



rp 



tfi) 



,(*) 



(20) 



where p« ,p^ e ' ,p^ p ' are the tetr ad(or LNRF) components 
of the four-momentu m (see e.g. IPineault fc Roederl Il977l ; 
IChandrasekharlll983h . 



2.4 Method of calculation 

We now turn to how to calculate the line profiles numeri- 
cally. A steady-state warped disc structure is first obtained 
by solving the dynamical equation numerically by setting 
dL/dt — in equation ©. The disc is divided into a num- 
ber of arbitrarily narrow rings, each such emitting ring being 
denoted by its radius n and weights uij, provided by an al- 
gorithm developed by Rybicki G. B. (|Press et al.lll992T ) for 
Gauss-Legendre integration. The main numerical procedures 
for computing the line profiles are as follows: 

(i) Using the numerical disc structure obtained first, the 
values of two Eulerian angles /3(R,t) and 7(i?,t) at an arbi- 
trary point are calculated by means of cubic spline interpo- 
lation. 

(ii) The relevant disc system parameters are specified: 
Rin, Rout, Ry, ^w, i, ">1 and 

(iii) For a given couple (n, g) of a ring, the two constants 
of motion A and q are determined if they exist. 

(iv) For each g, the integration over r of equation (|16l) 
can be replaced by 



F {u ) = J2 



qev a 



d(X,q) 



(21) 



where the emissivity e in equation (I2ip is calculated numer- 
ically by equation (|17[) . 

From the above formula, one can determine the line flux 
from the disc at an arbitrary frequency v . The observed line 
profile as a function of frequency v a is finally obtained in this 
way. 



3 RESULTS 

In the accretion disc model, the double-peaked emission 
lines are radiated from the disc region between around 
several hundred gravitational radii to more than 2000r g ; 
here, r g is the gravitational radius, and the widths of the 
double-peaked lines range from s everal thousand to nearly 
40,000 kms -1 (|Wang et al.ll2005l ). In our model, all of the 
parameters of the warped disc are set to be free. The warp 
radius i? w is a key parameter for descr ibing the disc shape. 
To det ermine the value of R w , either th ej shakura fc Sunvaevl 
(| 19731 ) or the lCollin fc Dumondl (|l990T ) model is used. As far 
as the outer region of the disc is concerned, the self-gravity 
of the disc must be taken into account at large radii. Be- 
yond a certain radius, the disc becomes gravitationally un- 
stable and can be disrupted. The typical range of unstable 
regio n is in the order of 10 3 to a few 10 4 Schwarzschild ra- 
dius (|Collin fc Hurdll999D . Thus, the reasonable range of 
the disc is from 100r 9 to 2000r g . In our calculation, the fre- 
quency ranges from 4.32 to 4.78 in units of 10 14 Hz, and 180 
bins are used. Considering the broadening arising from elec- 
tron scattering or turbulence, all our results are smoothed 
by convolution with a 3a Gaussian. 
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In (R/R») 

Figure 1. The distribution of f3 and 7 as a function of radius. 
The twist angle 7 is measured in units of 2-tt, and the tilt angle 
/3 is shown in radians. The solid, long-dashed, short-dashed and 
dotted lines show the results for 77 = 0, 3/4, 1 and 2, respectively. 



3.1 The numerical steady-state warped disc 
solutions 

We numerically solved equation ([5]), setting all the time 
differentials to zero, in order to find a steady-state disc 
structure. As a boundary condition, we set ft — ir/6 and 
7 = at the outer boundary. The calculation region is 
—4.6 < x < 2.3, where x = ln(R/R w ), or, equivalently, 
Rin = -R w e" 4 ' 6 « O.OLRw, and i? ou t = i? w e 2 ' 3 « 10i? w . The 
calculated distributions of f3 and 7 as a function of radius 
are plotted in Fig. [1] 

3.2 Line profiles from twisted warping discs 

We use the numerical code developed by IWu et alj (|2008h 
to evaluate the line profiles for the numerical disc struc- 
ture obtained. We first calculated the line profiles of a sim- 
plified model with no relativistic effects included in equa- 
tions of propagation of photons and compared them with 
the relativistic treatment. The parameters of the two Eule- 
rian angles an d the disc range are identical to those found by 
iBachevI (|l999l ) , but we used a ring- like primary source and a 
linear response for the reprocessing to the ionizing flux. The 
primary source radius is set to R r — 1.5r g . We reproduced 
the results obtained bv lBachevI l|l999T ) on a qualitative level. 
A comparison between the H/3 line profiles for a simplified 
model with no relativistic effects(solid line) and those with 
the ralativistic treatment (dotted line) is show in Fig. [2] The 
influence of the gravitational lensing, which concentrates the 
radiation from the innermost region towards the equatorial 
plane on the line profle, is significant. 

The dependence of the Ha line profiles on the pri- 
mary source radius Rr is ploted in Fig. |3] The radii 
with respect to the different lines are: R T = 1.5r g (solid 
line), 2. 0r g (dotted line), 3. 0r g (short-dashed line), 6. 0r g (dot- 
dashed line), 10. 0r g (long-dashed line). The other parame- 
ters are Ri n = 150r g , R out = 1200r g , i? w = 600r g , 97 = 1, 
1/2/vi = 10 and i = 30°. The azimuthal viewing angle along 
the XY-plane varies from 0° to 300° in steps of 60° from the 
top left-hand to the bottom right-hand panel. The vertical 
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Figure 2. Comparison of the H/3 line profiles between a simpli- 
fied model with no relativistic effects (solid line) and one with the 
ralativistic treatment (dotted line). The parameters of two Eule- 
rian angles an d the disc range are identical to those presented by 
Bachevl l|l999h . The azimuthal viewing angle along the XY-plane 
varies from 0° to 300° in steps of 60° from the top left-hand 
panel to the bottom right-hand panel. The vertical dashed line 
corresponds to the wavelength position of H/3 in the rest frame. 







n 


I ; 


- i \ 

- 1 } 




1 




■ /■I 


1 


1 




(trr 
1 


a\ 
1 








Q 

n 
\ 


rJ \ 





"wavelength 



Figure 3. Comparison of the Ho line profiles with different 
primary source radii R s : R s = 1.5r g (solid line), 2.0r g (dottcd 
line), 3.0r g (short dashed line), 6.0r g (dot-dashed line), 10.0r g (long 
dashed line). The other parameters are R ln = 150r g , i? ou t = 
1200r g , R w = 600r g , r) = 1, u 2 /vi = 10 and i = 30°. The az- 
imuthal viewing angle along the XY-plane varies from 0° to 300° 
in steps of 60° from top left-hand panel to bottom right-hand 
panel. The vertical dashed line corresponds to the wavelength 
position of Ho in the rest frame. 



dashed line corresponds to the wavelength position of Ha 
in the rest frame. From this figure, it can be seen that the 
line profiles are not sensitive to a change in the radius of the 
primary source in the range from 3.0r g to 10.0r g , but that 
they have a significant dependence when it approaches the 
innermost stable orbit. 

The images of the illuminated area and the Ha line 
profiles computed by our code for a steady-state twisted 
warped disc are shown in Fig. [4] The disc zone is from _Ri n = 
150r g to i?out = 1200r g , the warp radius is set to i? w = 
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Figure 4. Images (upper panels) of the illuminated area of the 
disc and the Ho line profiles (lower panels) computed by our code 
for the twisted warped disc case for i = 30° ; the disc zone is from 
R ln = 150r g to -Rout = 1200r g . The Bardeen-Petterson radius 
is set to i? w = 600r g , rj = 3/4 and vi/vx = 10. The azimuthal 
viewing angle along the XY-plane varies from 0° to 330° in steps 
of 30° from the top left-hand panel to the bottom right-hand 
panel. 
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Figure 6. Comparison of the Ho line profiles with various val- 
ues of the Bardcen-Pcttcrson radius or viewing angle. Upper 
two rows: the parameters are i = 30°, i?; n = 150r g ,i? O ut = 
2000r g , rj = 3/4, v 2 /ui = 10 and R„ = 1000r g (solid line), 
iiw = 500r g (dot-dashed line). Lower two rows: the parameters 
are R in = 150r g , R out = 1500r g , R w = 800r g , r) = 1, u 2 /vi = 10 
and i = 15° (solid line), i = 30° (dot-dashed line). The azimuthal 
angle along the XY-plane varies from 0° to 330° in steps of 30° 
from the top left-hand to the bottom right-hand panel. 
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Figure 5. Comparison of the Ho line profiles for a warped disc 
with various values of the power-law index of the viscosity: rj = 
0, 3/4, 1, 2 (from top to bottom rows). The other parameters are 
i = 30°, R in = 150r g , R out = 1200r g and R v = 600r g . The 
azimuthal viewing angle along the XY-plane varies from 0° to 
300° in steps of 60° (from the left- to the right-hand side). 

600r g , and the other parameters are i = 30°, ?7 = 3/4 and 
vilv\ = 10. The azimuthal viewing angle along the XY- 
plane varies from 0° to 330° in steps of 30° from the top left- 
hand to the bottom right-hand panel. The image contains 
400 x 360 pixels. From the images, we can see that, moving 
inwards from the outer radius of the disc, the warp is twisted 
by an angle of ~ ir before being flattened efficiently into the 
aligned plane, as is also shown in Fig. [T] 

The influence of the power-law index of the viscosity 
on the line profiles for r\ — 0, 3/4, 1, 2 (from top to bottom 
panels) is shown in Fig. [5] The disc zone is from R ln — 
150r g to i?out = 1200r g ; the other parameters are i = 30°, 



-R w = 600r g .The azimuthal viewing angle along the XY- 
plane varies from 0° to 300° in steps of 60° (from the left- to 
the right-hand side). The general behaviour of the line profile 
is similar. Fig.[6]shows the dependence of the profiles on the 
warp radius and viewing angle. The upper panels are for 
i = 30°, R m = 150r g , i? out = 2000r g , rj = 3/4, vs/vi = 10 
and i? w = 1000r g (solid line), -R w = 500r g (dot-dashed line). 
The lower panels show the profile dependence on the viewing 
angles for i = 15° (solid line) and i = 30° (dot-dashed 
line). The other parameters are R[ n — 150r g , -Rout = 1500r g , 
-Rw = 800r g , rj — 1 and va/vi = 10. 

Generally speaking, the profiles of the line emission 
from the warped disc induced by the Bardeen-Petterson ef- 
fect are nonsymmetrical and red- or blue-frequency-shifted 
in most cases for which the primary source radius is large. 
When the primary source radius approaches the innermost 
stable orbit, double-peaked or triplet-peaked line profiles 
are present in most cases. Thus, a warped disc induced 
by the Bardeen-Petterson effect provides a mechanism for 
producing the double- or triplet-peaked line profiles that 
were observed in the Sloan Digital Sky Survey. The triplet- 
peaked line profile, such as for SDSS J084 205. 57+075925 . 5 
and SDSS J232721.96+152437.3 (see e.g. IWu et ai]l2008h , 
may be a "signature" of a warped disc. Future monitoring 
of line profile variability in these triplet sources will provide 
a critical test of the warped disc model. 



4 SUMMARY 

We computed the Balmer emission-line profiles that result 
from the reprocessing of emission using a numerical model 
of a warped disc induced by the Bardeen-Petterson effect, 
including all relativistic effects for radiation propagation, 
and subject to various shadowing effects associated with 
disc warping. We numerically solved the disc structure for 
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a steady-state disc shape. To simplify the analysis, we took 
both viscosities to have the same power-law form, and thus 
the ratio vijv\ is independent of radius. For simplicity, we 
assumed that the disc is illuminated by a ring-like central 
source, which is a rough approximation, that the line emis- 
sivity is proportional to the continuum light intercepted by 
the accretion disc, and that line emission is isotropic. Our 
conclusions are as follows. 

(i) From the numerical results, we can see that, moving 
inwards from the outer radius of the disc, the warp is twisted 
by an angle of ~ 7r before being flattened efficiently into the 
equatorial plane. 

(ii) For a given warped disc, there are two angles that 
determine the observed line profile, the inclination angle 
and the azimuthal viewing angle resulting from the non- 
axisymmetry of the disc warping. 

(iii) For less twisted warped disc induced by the Bardeen- 
Petterson effect, the asymmetrical and frequency-shifted 
single-peaked line profiles are produced in most cases when 
the primary source radius is large(from 3r 9 to 10r g ). The 
double- or triplet-peaked line profiles presented in most 
cases occur when the primary source radius approaches the 
marginally stable orbit of an extreme Kerr black hole. 

The line prof iles depend stro ngly on the twisting struc- 
ture of the disc (|Wu et al.l 120081 ) as well as on the radius 
of the illuminating source. The Doppler effect and gravita- 
tional focusing concentrate the radiation from the innermost 
region towar ds the equatorial p l ane, w hich increases the re- 
turning flux l| Cunningham 19751. Il976h . This may be the rea- 
son why the illuminating source radius has a strong influence 
on the line profile. Knowledge of the twisting structure may 
allow the determination of some important characteristics of 
AGNs, such as the spin of the black hole and the viscosities 
of the disc. We assume that the ratio 1/2/ vi is independent of 
radius in the calculations, but the relation between v\ and 
V2 is very uncertain, and may influence the disc twisting 
structure. It is also important to emphasize that, although 
we have analyzed the Bardeen-Petterson effect in this work, 
the four mechanisms mentioned above are not mutually ex- 
clusive. This means that more than one mechanism might 
be operating in a given system. 
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